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Contributions from the reaction e + e~ — ► 7r + 7r~e + e~ to the pion form 
factor measurement via radiative return method are discussed basing on 
the results of a Monte Carlo generator (EKHARA). The generator con- 
tains contributions from the initial and final state emission of a e + e~ pair 
from e + e _ — > 7r + 7r~ production diagrams and the 7r + 7r~ pair production 
from space-like and time-like Bhabha diagrams. A detailed study is per- 
formed for the $- factory energy. Tests of the generation procedure are 
also presented. 

PACS numbers: 13.40.Ks,13.66.Bc 

1. Introduction 

The radiative return method is a powerful tool in the measurement of 
a(e + e~ — ► hadrons) and detailed studies of hadronic interactions [21 El II]- 
Very accurate knowledge of the hadronic cross section is essential for predic- 
tions of the hadronic contributions to cu, the anomalous magnetic moment 
of the muon, and to the running of the electromagnetic coupling from its 
value at low energy up to Mz (for recent reviews look [HI El El)- Due 
to a complicated experimental setup, the use of Monte Carlo (MC) event 
generators [HI [HI , which include various radiative corrections jlUj is indis- 
pensable. Some more extensive analysis of that subject can be found also in 
jllj . The most important hadronic mode, i.e. 7r + 7r~, was recently measured 
by KLOE ^2] by means of radiative return method. In this measurement 
only pions (+ missing momenta) in the final state were observed. For that 
particular measurement there is no difference between photon(s) and the 
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pair production and one should estimate the contribution of the process 
e + e~ — > 7r + 7r - e + e~ to the measured cross section. As suggested in ^3] the 
contribution from the e + e~ production to the e + e~ — > 7r + 7r~7 is sizable 
and comes mostly from the t-channel Bhabha-like diagrams. In this paper 
a Monte Carlo study is performed to test if this claim remains true for a 
realistic (KLOE) experimental setup. 

2. Monte Carlo simulation and its tests 

a) b) 




Fig. 1. Diagrams contributing to the process e + (pi)e~(p2) — > 

7r + (7r 1 )7r _ (7r2)e + (qi)e _ ((72): initial state electron pair emission (a), final 
state electron pair emission (b), pion pair emission from t-channel Bhabha process 
(c), 7*7* pion pair production (d) and pion pair emission from s-channel Bhabha 
process (e). 

In Fig^different types of diagrams contributing to the reaction e + e~ — > 
7r + 7r~e + e~ are shown schematically. In the present version of the Monte 
Carlo program we include all but the 7*7* pion pair production process 
diagrams (Fig^l), which were estimated to be negligible for DA<I>NE energy 
|14j . We use scalar QED to model the FSR e + e~ pair emission and the 
p dominance model for the 7*(/0*)7T7t coupling (see E3 and ^H] (from 
where the model of the pion form factor was implemented) for details). 
The contribution from initial state radiation (ISR) diagrams Fig. and 
final state radiation (FSR) diagrams Fig. ^> was discussed in details in jUJ . 
Because of a heavy pair emission, the contribution from diagrams presented 
schematically in Fig. ^ is completely negligible for any event selection used 
in the analysis presented below. 
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We use multi-channel variance reduction method to improve efficiency 
of the generator and the generation is split into four channels, where two of 
them absorb peaks present in t-channel diagrams and other two take care 
of the s-channel peaks. All details will be given in a separate publication 

na. 

To check correctness of the program code, we have performed a number 
of tests of the new part of program. Gauge invariance of the sum of the 
amplitudes was checked analytically for set of diagrams from Fig. ^,e. To 
cross check the helicity amplitudes, used in the program to calculate square 
of the matrix element, we have used also the standard trace method for inde- 
pendent calculation. Both results, summed over polarizations of initial and 
final leptons, were compared numerically scanning the physical phase space. 
The biggest relative difference between the two results, which was found, 
was at the level of 10 -11 for diagrams in Fig. ^ and 10 -23 for diagrams in 
Fig. The computer code has been written in quadruple precision not to 
lose accuracy as sever numerical cancellations occur even when using helic- 
ity amplitudes. The phase space volume as a function of the square of the 
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Fig. 2. Relative difference between phase space calculated by MC and by gauss 
method. The errors come from MC integration (left). Relative difference between 
MC result and average for phase space x simple \M\ 2 (right). 

invariant mass of the pion system calculated by the Monte Carlo program 
was compared with a one dimensional Gauss integration. The result is pre- 
sented in Fig. The relative difference stays well within errors, which are 
at most at the level of a few per mil. Because of the complexity of the gen- 
eration we performed an additional test where a very simple matrix element 
was introduced. The matrix element was chosen in a way that only a four 
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Table 1 . Ratios Rt+int and R in % (see text for definitions) : BN — results based 
on formulae from BNK — results based on formulae from JJj and ^H], MC 
— Monte Carlo results 



dimensional numerical integration was necessary. The matrix element was 
integrated by the Monte Carlo program with full phase space generation and 
compared with the results obtained with 4-dimensional Gauss integration, 
where the integration region was divided into (2 • 10 7 — 6 • 10 7 ) rectangular 
solids (or more complicated shapes at the border of integration region). In 
each rectangular solid the 8-point Gauss procedure was used recursively. 
The error of the Gauss method was estimated as the maximal difference 
of the results with different subdivisions and the central value was taken 
as a mean of the smallest and the biggest obtained values. The relative 
difference between MC result and the Gauss method is presented in Fig|2J)- 
Again the relative difference stays well within errors, which are at most at 
the level of a few per mil. 

A comparison with existing analytic results |17[ I18| was also performed. 
In the Tabled the ratio (denoted by Rt+int) of differential cross sections for 
t-channel (space-like Bhabha with radiated pion pair) plus the interference 
between t-channel and s-channel ( only the ISR pion pair radiated from 
e + e~ — > e + e~ is included) to the differential cross section of the reaction 
e + e~ — > 7r + 7r~7 (Born: ISR+FSR) is presented as a function of the invariant 
mass of the two-pion system (Q 2 ). No cuts are imposed as the analytic 
formulae exist only for that case. A discrepancy at the level of about 10 
% between the Monte Carlo and the analytic results can be observed in 
the region were the pair contribution is sizable and it is much larger in the 
region were the pair contribution is at the level of 1 per mil or lower. For 
s-channel diagrams we have found previously jH] an agreement with analytic 
results of ^S] at the level of 1 per mil. We shall observe, that even if we 
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take the relative difference as the error of the MC code, the accuracy of 
the code is good enough for an estimation of the contribution of the pair 
production to the process with the photon(s) radiation. To estimate the 
complete contribution for the (unrealistic) situation with no cuts imposed, 
the ratio (denoted by R) of the t- and the s-channel pair production to 
the differential cross section of e + e~ — > 7r + 7r~7 (Born: ISR+FSR) is also 
presented in Table^both for the MC based and the analytic results [171 118] . 



3. An event selection dependence of the pion pair production 
contribution to the radiative return measurement 
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Fig. 3. (a) Ratio of differential cross sections of t-channel and s-channel of the 
reaction e + e~ — > n + ir~ e + e~ (b) Ratio of differential cross sections of the reactions 
e + e~ — ► 7r + 7r~e + e~ (pair) and e + e~ — > 7r + 7r~7 (ph) 

All results from Monte Carlo simulation, presented in this section, are 
for DA<&NE energy 1.02 GeV. As observed already in the paper ^3] t- 
channel diagrams contribute significantly in the small Q 2 region if no cuts 
are imposed. That is shown in Fig. |31 where Fig. presents the ratio of 
t-channel to s-channel contributions to the differential cross section of the 
reaction e + e~ — > ir + Tr~e + e~ and Fig.|3p presents the ratio of the differential 
cross sections of the reactions e + e~ — > 7r + 7r~e + e - and e + e~ — > tt + tt^j. 
The t-channel diagrams contribution, if no cuts are imposed, is up to 12 
times bigger than the s-channel contribution in the small Q 2 region and it 
is negligible for large Q 2 values. 

The contribution from the pair emission to the photon emission depends 
however a lot on the event selection, even if the t-channel contribution still 
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Fig. 4. Ratio of differential cross sections of the process e + e~ — > 7r + 7r~e + e _ (pair) 
and e + e~ — > 7r + 7r~7 (ph) with two sets of cuts imposed (a) KLOE cuts (b) only 
pion angular cuts are imposed 
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Fig. 5. Charge asymmetry ^(e^>9o°) + N(fl^<9o°j : ( a ) for e+e ~* 7r+7r (7 + 
e + e~) (the solid line shows the charge asymmetry for e + e~ — > 7r + 7r~7 only) and 
(b) for e + e~ — > 7r + 7r~e + e~ (no cuts imposed) 



dominates at low Q 2 . It is seen from Fig. 0J where for KLOE event selection 
(Fig. Efci) one gets only up to 1.2% contribution from the pair emission and 
up to 2.5% if one uses only the angular cuts on pions polar angles (Fig. HJa). 
The pair contribution at the p peak is about 0.3% for KLOE cuts and 0.5% 
when only the angular cuts are imposed. Even if the pair contribution 
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Fig. 6. Charge asymmetry ;v(e"+>9o°)+Ar(^+<9o°) : for ( a ) e+e 7r+7r (7 + 
e + e _ ) (the solid line shows the charge asymmetry for e + e~ — > 7r + 7r~7 only) and 
(b) for e + e~ — » 7r + 7r _ e + e _ (for KLOE event selection) 



remains small it has to be taken into account if the aimed experimental 
precision is of the order of 1% or better. 

Another important observable used to control the FSR contribution to 
the cross section of the reaction e + e~ — > 7r + 7r~ + photons is the charge 
asymmetry. Again the measured value contains the pair contributions if 
only pions are observed. However one can see in Figs. EH and EH that 
this contribution is completely negligible both in the case where no event 
selection is used (Fig. EH) and in the case when KLOE event selection is 
applied (Fig. EH)- That conclusion holds despite the fact that the asymmetry 
for pair production itself is sizable in both cases (Figs. Eb andEb)- 

As stated already in j^, ISR and t-channel of electron pairs can be 
treated likewise ISR of photons, what results in the change of the radiator 
function in the radiative return method. Another possibility is to treat the 
pair production as a background to the process e + e~ — > 7r + 7r~7(7). In that 
case one would like to reduce that background. The way to do that is to 
veto the outgoing electrons and/or positrons. However as shown in Fig. [7[ 
where the percentage of the electrons and positrons escaping the detection 
((0 e + < 20° or 6 e + > 160°) and (0 e - < 20° or 9 e - > 160°)) is plotted, it is 
clear that it is practically impossible for KLOE event selection (Fig. EH), as 
about 90% of electron-positron pairs escapes detection. For looser cuts one 
can reduce up to 50% of the pair background (Fig. da). 
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Fig. 7. Non-reducible pair production background 



4. Conclusions 

A detailed analysis of the contribution of the reaction e + e~ — > ir + ir~e + e~ 
to the e + e~ — > 7r + 7r _ cross section extraction via the radiative return 
method from the a(e + e~ — > 7r + 7r _ + missing momentum) measurement is 
presented basing on the Monte Carlo simulations with the EKHARA event 
generator. It is shown that t-channel diagrams contributions dominate at 
small Q 2 values independently on the event selection. The biggest possible 
contribution to the process e + e~ — > 7r + 7r~7 reaches about 8 % if no event 
selection is applied and it is reduced to 1.2 % for KLOE event selection. 
The reaction e + e~ — > 7r + 7r~e + e - does not contribute significantly to charge 
asymmetry measured in the reaction e + e~ — > 7r + 7r~ + missing momentum. 
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